Soliton Magnetization Dynamics in Spin-Orbit Coupled Bose-Einstein Condensates 
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Ring-trapped Bose-Einstein condensates subject to spin-orbit coupling support localized dark 
soliton excitations that show periodic density dynamics in real space. In addition to the density 
feature, solitons also carry a localized pseudo-spin magnetization that exhibits a rich and tunable 
dynamics. Analytic results for Rashba-type spin-orbit coupling and spin-invariant interactions pre- 
dict a conserved magnitude and precessional motion for the soliton magnetization that allows for 
the simulation of spin-related geometric phases recently seen in electronic transport measurements. 
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The recent realization of artificial light-induced gauge 
potentials for neutral atoms [1 has added a power- 
ful new instrument to the atomic-physics simulation 
toolkit [2]. In particular, possibilities to induce Zeeman- 
like and spin-orbit-type couplings in (pseudo-)spinor 
atom gases [3 render them ideal laboratories to investi- 
gate the intriguing interplay of spin dynamics and quan- 
tum confinement that has been the hallmark of semicon- 
ductor spintronics [H |5] . At the same time, the unique 
aspects of Bose-Einstein-condensed atom gases [6 asso- 
ciated, e.g., with their intrinsically nonlinear dynamics, 
promise to give rise to novel behavior under the influence 
of synthetic spin-orbit couplings [TWTS] . 

One of the special properties resulting from nonlinear- 
ity in Bose-Einstein condensates (BECs) is the existence 
of solitary- wave excitations [16]. Basic types of these 
are distinguished by the shape of their localized density 
feature: dark (gray) solitons are associated with a full 
(partial) depletion of a uniform condensate density in a 
finite region of space, whereas bright solitons are local- 
ized density waves on an empty background. A further 
characteristic associated with solitons is the phase gra- 
dient of the condensate order parameter centered at the 
position of the density feature. In multi-component sys- 
tems, the dynamics of soliton excitations is found to be 
enriched by the additional degrees of freedom p!7H2Q] . 

We have studied solitons in ring-trapped pseudo-spin- 
1/2 condensates with spin- invariant repulsive atom- atom 
interactions subject to a Rashba-type [211 122] spin-orbit 
coupling and find that they exhibit a third feature: a 
pseudo-magnetization vector with conserved magnitude 
and rich dynamics that unfolds in tandem with the soli- 
ton's periodic propagation in real space. Figure [T] shows 
an example and also illustrates the interesting fact that 
the magnetization directions at the beginning and the 
end of a full cycle of the soliton's motion are generally not 
parallel. The appearance of such a geometric phase [23] 
and the precessional time evolution of the solitonic mag- 
netization is reminiscent of the spin dynamics of electrons 
traversing a mesoscopic semiconductor ring 



In the following, we consider several soliton configura- 
tions and obtain analytical results for their density and 
magnetization profiles as well as the magnetization dy- 
namics associated with their motion. We start by intro- 
ducing the basic theoretical description of our system of 
interest. Using the basis of a spatially varying local spin 
frame [26] for the condensate spinor, the nonlinear Gross- 
Pitaevskii equation [6 for the spin-orbit-coupled ring 
BEG turns out to be of Manakov-type [17] , making it pos- 
sible to apply standard methods [I8l [20] to find solitary- 
wave solutions. Accounting for the presence of spin-orbit 
coupling adds an important twist: Spinors have to sat- 
isfy non-standard boundary conditions, which introduce 
background-density fiows in the local spin frame that 
contribute to the nontrivial magnetization dynamics ex- 
hibited by the moving solitons in the lab frame. 

We consider a two-component BEG trapped in the xy 
plane and confined to a ring of radius R. The atoms 
are assumed to be in the lowest quasi-onedimensional 




FIG. 1. Time evolution of a gray-bright soliton's magnetiza- 
tion in a ring-trapped EEC with Rashba spin-orbit coupling. 
During a full cycle of the soliton's motion on the ring, the 
magnetization vector follows a trajectory on the surface of a 
sphere. The magnetization vectors at the beginning and the 
end of a cycle (indicated by arrows) differ by an angle i9 that is 
related to a spin-related geometric phase. Soliton parameters 
(see text): Vs/c = 0.5, tan 77 = 2, ^ = 100, k, = —0.01. 



subband [29] and subject to a spin-orbit coupling of the 
familiar Rashba form [lU |22] aYi[ax{—idy) — cry{—idx)] 
as well as a spin-rotationally invariant contact inter- 
action. (Here crx,y,z are the spin- 1/2 Pauli matrices.) 
The energy functional of such a system [30] is given by 
E[^] = J dip "^^ {H — /i) ^, where (p is the azimuthal 
angle, fi the chemical potential, ^ = ('0|,'^^)^ the two- 
component (pseudo-spin- 1/2) spinor order parameter in 
the representation where the (z) direction perpendicular 
to the ring's plane is the spin-quantization axis, and 
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We use a± = {ax ± ^o-y)/2 to denote raising and lower- 
ing operators for spin-1/2 components, Eq = h^ /{2MR^) 
is the energy scale for quantum confinement of atoms 
with mass M in a ring of radius i?, E^g is the two-body 
contact-interaction strength, and tan 7^ = 2MRaYi/h^ is 
a dimensionless measure of the spin-orbit coupling. 

The effect of Rashba spin-orbit coupling in a ring ge- 
ometry can be elucidated by performing a suitable SU(2) 
transformation. Defining ^ =Ux and l-L\oc = U~^1-LU^ 

with U = Q-^^^zl1 ^ir)cjyl2 ^i^cj,l{2co^r)) ^^ ^^^ 



u 



loc 



^n 



-dl 



(tan 7^)^ 



9_ t 



(2) 



The transformation IA~^ amounts to a (/^-dependent rota- 
tion of the pseudo-spin quantization axis [26 , followed by 
a spin-dependent gauge transformation. We will refer to 
the original representation where the spin-quantization 
axis coincides with the axis of the ring as the lah frame ^ 
whereas the representation in which the Hamiltonian of 
the system is diagonal in pseudo-spin space [i.e., given 
by Hioc of Eq. ([2|] will be the local spin frame [26 . Note 
that the spinors ^ in the lab frame are periodic func- 
tions of (/?, whereas the spinors x — (x+^X-)^ from the 
local spin frame have to satisfy the boundary conditions 
X±(^) = X±(^ + 27r)e^*'^ with a spin dependent phase 
twist originating from the spin-orbit coupling, where 
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Knowledge of the local-spin-frame spinors enables the 
calculation of expectation values for any observables ac- 
cessible to measurement in the lab frame. The total 
density n = IV^^P + IV^^P = |x+P + IX-P is obviously 
the same irrespective of which representation is chosen in 
spin space. The pseudo-spin- 1/2 projections in the lab 
frame correspond to definite atomic states, hence their 
density profiles n^{i) = ^^([1 + {—)az]/2)'^ are of in- 
terest. In addition, we will consider the magnetization- 
density vector s = "^^a"^ in the lab frame, with a = 
{ax^ (Ty,cFz) being the vector of Pauli matrices. 



We analyze the properties of localized excitations in 
spin-orbit-coupled ring-trapped BEC based on the time- 
dependent Gross-Pitaevskii equation [6 5E[x\/5x*(j = 
ihdxa/dt' After rescaling to use the dimensionless time 
variable r = tEo/h^ it has the form 
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for the two components of the spinor x = (x+^X-)"^? 
where no = [/i + (tan77)^/4]/(^£^o) is the uniform (back- 
ground) density consistent with the chemical potential 
(1. While the spin-orbit coupling has formally disap- 
peared from the nonlinear equation Q, it is still im- 
plicitly present via the boundary conditions that the in- 
dividual components X±{^->'^) must satisfy. 

We have obtained several soliton solutions of Eqs. Q 
using established techniques [T7l[T8l[20] and implemented 
the appropriate boundary conditions. Before giving fur- 
ther details, we like to summarize a few general features. 
The soliton spinors in the local-spin-frame representation 
turn out to be of the form 



^(s) ^ ^^^){ip - y^r) e*^b^^/2-iVb<T^/4 
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where T^ (0 are complex amplitude functions encoding 
the specific soliton-like density features, v^ is the propa- 
gation speed of the soliton, and v\^cr are background flow 
velocities of the individual spinor components that are 
necessary to implement the boundary conditions arising 
due to the presence of spin-orbit coupling. The density 
and magnetization density s^®^ exhibit spatially lo- 






calized features. Subtracting s*^^^ from the magnetiza- 

(s) 

tion density s^ of the condensate background yields the 
magnetization density that is associated with the soli- 
ton excitation only. Its integral S^^-* = J dcp [s\^^ — s*^^^] 
is the vector of total soliton magnetization, which is 
an additional property of localized excitations in multi- 
component BECs. For soliton solutions of the form ([5|, 
S^^-' has constant magnitude. Its temporal evolution is 
most conveniently described by a set of four angles as 
defined in Fig. [2|c) . While the tilt angles (3 and (3' are 
time-independent, the angles a and a' vary linearly in 
time, signifying the precession of S^^^ around tilted z^ 
axis with the universal result 



13 = V, 



VsT + TT. 
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The z^ axis is tilted by the angle r] characterizing the 
spin-orbit coupling and it rotates around the z axis with 
the same angular velocity v^ that characterizes the soli- 
ton propagation. The second tilt angle /3' is found to de- 
pend only on the soliton profile T± ((/:?), while the preces- 
sion frequency da' /dr has complicated dependences on 
the parameters of the soliton solutions. Figure |2] shows 
exemplary magnetization dynamics for gray-bright and 




FIG. 2. Magnetization dynamics of a gray-bright soliton [panel (a)] and a gray- gray soliton with zero background magnetization 
in the local spin frame [panel (b)]. Parameters used are g = 100 (100), tan 77 = 0.2 (0.5), Vs/c = 0.5 (0.2), and k = —0.5 for 
the gray-bright (gray-gray) case, (c) Angles used to describe the two-step precessional motion of S. The angle a is measured 
with respect to an x' axis that is perpendicular to both the z and z' axes. 



gray-gray solitons. Interestingly, we find that the mag- 
netization vector is usually not parallel to its initial di- 
rection after the soliton has completed a full cycle of its 
motion around the ring as, e.g. seen in figure [l] The an- 
gle 1^ between the magnetization directions at the start 
and the end of a cycle turns out to be finite only as a 
consequence of spin-orbit coupling, as it depends promi- 
nently on the phase A given in Eq. ([3| that also governs 
spin-dependent interference in mesoscopic ring conduc- 
tors [27 . 

In order to find explicit soliton solutions, we introduce 
^ = (p — ur, where ix is a velocity parameter, and initially 
look for solutions of the form Xcr(v^,T) = \Av(0^*' 
This allows us to rewrite Eq. Q in the form 
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where n = U-^ -\- U-. Single component solutions for 
a = a are easily found by integration of Eqs. ^ to yield 
Xa(0 ^ ^(0 ^^^ X-a = 0, with the well-known dark 
soliton solution on the infinite line 6 
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1- +7ntanh 7^ -— 
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Here 7^ = 1 - u^/c^, c^ = 2gno, l/^l = gno/2. The soli- 
ton profile ([8| is appropriate for sufficiently strong non- 
linearity, where (,d/1u <^ 27r pTj. However, T(^) does 
not satisfy the proper boundary condition since it has a 
phase step A^ = — 2arccos(ii/c). To compensate for the 
phase step and ensure the correct phase shift associated 
with the gauge transformation U^ we perform a Galilean 
transformation on ([8|, which yields 



xT{^.r) 



T(^ 
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Here v\) = —{A6 + aA)/7r is the background velocity 
imposed by the boundary condition. Thus the single- 
component soliton solution is of the form (l5|, with T|^ = 
T and T^^ = 0, and propagation speed Vs = u -\- vi^. 



A straightforward calculation yields s^^ = a\T{(p — 
Vsr)\'^{— sin T] cos (p, — sin T] simp ^ cos T])^ for the magne- 
tization density of the single-component soliton solution. 
In essence, the density depletion at the soliton's posi- 
tion gives rise to a reduction of the magnetization den- 
sity s^^ = crno(— sinr^ cos(/?, — sinr^ sin(/:?,cos7^)^ asso- 
ciated with the background. Thus s^^ — s^^ consti- 
tutes the magnetization density associated with the soli- 
ton itself, as it is the change in the background mag- 
netization density due to the presence of the localized 
excitation. For the single-component soliton, this corre- 
sponds to a peak in magnetization density at the soli- 
ton's position. The total magnetization vector is ob- 
tained by integrating that peak in real space, which yields 
S = cr(— sin?] cos(^'s'?"), — sinr^ sin('L'sT), cos?])-^. This 
magnetization vector is precessing in a perfectly synchro- 
nized fashion with the soliton's motion around the ring 
[cf. Fig.jSJc) with Eq. ^ and P' = 0], i.e., i9^^ = 0. 

We now consider a solution of Eqs. ^ that is a gray- 
bright (GB) soliton in the local spin frame. We assume 
that the densities approach constant values n+ -^ no+ 
(gray part) and n_ ^ (bright part) far away from 
the soliton's position. To decouple Eq. (7b), we use the 
ansatz [20 n_ = A>:(n+ — no+) with — 1 < /^ < 0. We 
apply a Galilean boost to both components to match the 
phase of the gray part only, hence they are of the form 
^ with Tf^{^) given by T(0 from Eq. ^ but with 
rescaled c^ = 2^no+(l + tz)^ 1/^^ — 5'^o+(l + ^)/2, and 



Te^(0 = V^^^ 
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cosh(7^^/^D) 
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Furthermore, v\^- = vi^-^ and v^ = u -\- vi^-^. Figure [3] 
shows the density profiles [panel (a)] and magnetization- 
density profile [panel (c)] associated with a GB soliton. 

The vector S of total magnetization for a GB soli- 
ton precesses concomitantly with the soliton's motion; cf. 
Fig.|2];c) withEq. ^ and tan /3' = {^/^U7T)/[{l-l^)c-fu], 
a' = —Vsr{l -\- A/tt). Figures flla) and[2|a) show exam- 
ples of possible time evolutions of the GB-soliton magne- 
tization. The magnitude of the magnetization vector is 



1.0 
0.8 
0.6 
0.4 
0.2 



:La)__ 


"f""-- 




1 

/ \1 ■ 


^ """"' '!■/ i 



0.2 0.4 0.6 0.8 



1.0 




1.0 
0.5 
0.0 



-0.5 



(c) 



^ I 

!\ P, 

\\ ' 
W 



1.0 



0.2 0.4 0.6 0.8 1.0 



(d) 



0.5 

0.0 

-0.5 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 3. Lab- frame spin densities [(a),(b)] and z compo- 
nent of the magnetization density [(c), (d)] for a stationary 
gray-bright soUton [(a),(c)] and a stationary gray-gray soUton 
[(b), (d)]. Panels (a) and (b) show the total density (solid yel- 
low curve) and individual-spin (dashed blue = t? dotted red 
= \) densities normalized to the background- density value. In 
(c) and (d), the total density is plotted again for reference as 
the solid yellow curve, together with the magnetization pro- 
files Sz (dashed blue curve), and the magnetization density 



» 



Sz associated with the soliton excitations only (dotted 



red curve) . Parameters are tan 77 : 
gray-bright soliton) n = —0.5. 



1.0, ^= 100, and (for the 



found to be 5* = 2(1 — n)(^D"fu/ ^^^ P' - For a GB soliton, 
the magnetization vector turns out to be not aligned with 
its initial direction after completion of a full cycle of its 
motion around the ring. A straightforward calculation 
yields sin('i^^^/2) = sin/3' sin^. As j3' is a known func- 
tion of the soliton parameters, a measurement of 1?^^ will 
yield the spin-related geometric phase A. 

The solutions representing gray-gray (GG) solitons in 
the local spin frame are obtained by Hirota's method [18 . 
The spinor components are of the form ([5| with 
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Here 

^2 _ 
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gnocrf^- The back-ground flows 
are given by I'bcr = — (A^^r + (tA)/7t^ where AOo- = 
— 2arccos(iicr/ccr). The independent parameters charac- 
terizing a GG soliton are the ratio no+/no- (or, equiv- 
alently, the background magnetization in the local spin 
frame) and the speed v^ of the soliton. All other param- 
eters can be found by solving transcendental equations 



given just after Eq. (11). For simplicity, we consider the 



case of a GG soliton with zero background magnetiza- 
tion in the local spin frame (i.e., no+ = no- = no/2). 
Figure [3] shows results for spinor-density [panel (b)] and 
magnetization-density [panel (d)] profiles. 

The time evolution of the magnetization vector asso- 
ciated with a moving GG soliton is characterized by the 



angles defined in Fig. [2|c) with Eq. ([6|, /3' = 7r/2, and 
a^ = ur -\- 7r/2, where u = — ('^b+ — '^b-)'^s/2 + {v^_^ — 
vl_)/4:—{l^A/7r)vs. Figurel2lb) illustrates this dynam- 
ics which, for small 7^, corresponds to a slow rotation of 
the magnetization vector in the ring's plane with super- 
imposed fast small-amplitude oscillations in the normal 
direction. As in the case of the GB soliton, the magneti- 
zation vector does not evolve back to its initial direction 
after a period of the soliton's ring revolution. The angle 
between magnetizations at the start and the end of the 
cycle is found to be ^^^ = 27t\uj\/vs -^ 2^ for ^ <C 1. 
Again, the dependence of "d^^ on A enables determina- 
tion of the latter by measuring the former. 

In conclusion, we have investigated the properties 
of soliton excitations in ring-trapped spin-orbit-coupled 
BECs. We find that a magnetization degree of freedom is 
generally associated with a soliton, and that the magne- 
tization vector precesses around an axis that is rotating 
synchronously with the soliton's orbital motion around 
the ring. The magnetization direction at the end of a 
cycle of revolution does not coincide with the initial di- 
rection for the gray-bright and gray-gray cases, making it 
possible to measure a spin-orbit-related geometric phase. 
Our work opens up new avenues for the realization and 
manipulation of magnetic soliton excitations in BECs. 
It also creates the opportunity to study spin-dependent 
interference and scattering effects that, until now, were 
only accessible in semiconductor nanostructures. 
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